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Abstract 

We investigate the wave properties of the Kerr black hole with 
isothermal plasma using 3+1 ADM formalism. The corresponding 
Fourier analyzed perturbed GRMHD equations are used to obtain the 
dispersion relations. These relations lead to the real values of the 
components of wave vector k which are used to evaluate the quanti- 
ties like phase and group velocities etc. These have been discussed 
graphically in the neighborhood of the pair production region. The 
results obtained verify the conclusion of Mackay et al. according to 
which rotation of a black hole is required for negative phase velocity 
propagation. 

Keywords: 3+1 formalism, GRMHD, Kerr planar analogue, isothermal 
plasma. 

1 Introduction 



The study of gravitomagnetic waves in rotating black hole is important be- 
cause the existence of black holes can be ultimately verified with the help 
of infalling plasma radiation and super-radiance of gravitomagnetic waves. 
Information from the magnetosphere can be transmitted from one region to 
another only by means of propagation allowed by plasma state. The gravita- 
tional field and the wind bring perturbations to the external fluid dynamics. 
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The response of black holes to the external perturbations can be explored 
with the help of wave scattering method. 

General Relativity is the theory of four- dimensional spacetime but we 
experience a three-dimensional space evolved in time. It is easier to split the 
spacetime into three-dimensional space and one-dimensional time to develop 
a better understanding of the physical phenomenon. The split we usually use 
in understanding the general relativistic physics of black holes and plasmas is 
3+1 ADM split introduced by Arnowitt et al. pQ. This split in the formala- 
tion of general relativity is particularly appropriate for applications to the 
black hole theory as described by Thorne et al. [2]-[lj. Using this formalism, 
the wave propagation theory in the Friedmann universe was investigated by 
Holcomb and Tajima [5], Holcomb [B] and Dettmann et al. [7]. Komissarov 
|8j discussed the famous Blandford-Znajek solution. 

Blandford and Znajek [9J found a process which describes the extraction 
of rotational energy in the form of Poynting flux. The black hole with a 
force free magnetosphere behaves as a battery with internal resistivity in a 
circuit made by poloidal current. This current system is considered to be 
equivalent to incoming and outgoing waves. The incoming waves transport 
energy in a direction opposite to the Poynting flux. Penrose [10] was the 
pioneer who gave the idea of extraction of energy from the rotating black 
hole by a specific process called Penrose process. In the wavelength analogue 
of Penrose process [UJ an incoming wave with positive energy splits up into a 
transmitted wave with negative energy and a refracted wave with enhanced 
positive energy. The negative energy wave propagates into the black hole 
equivalent to a positive Poynting flux coming out of the horizon [12j-[13j. 

The key features of the Kerr black hole were beautifully summarized by 
Miiller [13] who investigated the accretion physics in the plasma regime of the 
general relativistic magnetohydrodynamics (GRMHD). Punsley et al. [T3] 
considered the black hole magnetohydrodynamics in a broader sense. Musil 
and Karas [16J observed the evolution of disturbances originated in outer 
parts of the accretion disk and developed a numerical scheme to show the 
transmission and reflection of waves. Koide et al. [T7] modeled the GRMHD 
behavior of plasma flowing into rotating black hole in a magnetic field. They 
showed (numerical simulations) that energy of the spinning black hole can 
be extracted magnetically. Zhang [18]- [19] formulated the black hole theory 
for stationary symmetric GRMHD with its applications in Kerr geometry. 
He discussed wave modes for cold plasma with specific interface conditions. 
Buzzi et al. [2D]- [21] provided a linearized treatment of transverse and lon- 
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gitudinal plasma waves in general relativistic two component plasma (3+1 
ADM formalism) propagating in radial direction close to the Schwarzschild 
horizon. 

Mackey et al. [22] gave the idea that negative phase velocity plane wave 
propagates in the ergosphere of a rotating black hole. They verify that the 
rotation of a black hole is required for negative phase velocity propagation 
which is a characteristic of Veselago medium. This medium was hypotheti- 
cally mentioned by Veselago [23] and later formed experimentally [23] as a 
material (called metamaterial or left-handed material). Much work has been 
carried out to investigate the characteristics of this medium [25]. Woodley 
and Mojahedi [2H] showed (using full wave simulations and analytical tech- 
niques) that in left-handed materials, the group velocity can be either positive 
(backwards wave propagation) or negative. Sharif and Umber [27]- [28] inves- 
tigated some properties of plasma waves by investigating real wave numbers. 
The analysis has been done for the cold as well as isothermal plasmas living 
in the neighborhood of the event horizon by using Rindler approximaton of 
the Schwarzschild spacetime. In a recent paper, the same authors [2H] have 
found some interesting properties of cold plasma waves using perturbation 
wave analysis to the GRMHD equations in the vicinity of the Kerr black 
hole. They have also discussed the existence of Veselago medium near the 
pair production region. 

This paper has been extended to investigate the wave properties for the 
isothermal plasma. We have focussed this work to investigate the following 
three main objectives: 

1. The behavior of gravitomagnetic waves under the influence of grav- 
ity and magnetospheric wind is analysed. This helps us to detect the 
response of the black hole magnetospheric plasma oscillations to grav- 
itomagnetic perturbations near the pair production region. The pair 
production region lies near the event horizon of the black hole. 

2. The existence of Veselago medium in the black hole regime is checked. 

3. The negative phase velocity propagation regions are investigated and 
compare the results with those obtained by Mackay et al. |22j . 

To this end, we derive the GRMHD equations in 3+1 formalism using the 
isothermal equation of state. The component form of the equations for spe- 
cific background assumptions is obtained by using perturbations. We con- 
sider the perturbed quantities as plane harmonic waves produced by gravity 
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and wind due to black hole rotation. The Fourier analysis method for waves 
is applied and dispersion relations are derived. These relations lead to the 
x-component of the wave vector from which the relevant quantities are in- 
vestigated to analyze the wave properties near the pair production region. 

The paper is organized as follows. The next section is oriented with 
the description of the Kerr analogue spacetime and mathematical laws in 
3+1 formalism for this model. Section 3 is devoted to the assumptions 
corresponding to the background flow. In section 4, the GRMHD equations 
alongwith their Fourier analyzed perturbed form for the isothermal equation 
of state of plasma are given. Section 5 provides the solutions of dispersion 
relations. In the last section, we shall discuss the results. 

2 Mathematical Framework 

This section contains the line element for a general spacetime model. The 
electrodynamics corresponding to Kerr planar analogue in 3+1 formalism is 
also considered. 

2.1 Description of Model Spacetime 

The line element of the spacetime in 3+1 formalism can be written as 

ds 2 = -a 2 dt 2 + 7 i i ((2x i + (3 i dt){dx j + (3 j dt), (2.1) 

where lapse function (a), shift vector (/3) and spatial metric (7) are functions 
of time and space coordinates. 

We consider the planar analogue of Kerr spacetime [19], i.e., 

ds 2 = -dt 2 + {dx + (3(z)dt) 2 + dy 2 + dz 2 . (2.2) 

Here z, x, y and t correspond to Kerr's radial r, axial 0, poloidal 6 and 
time t coordinates. The Kerr metric depends non-trivially on both r and 6, 
whereas this model metric depends on z only. The lapse function a is taken 
to be unity to avoid the effects of horizon and redshift. The value of the shift 
function (3 (analogue to the Kerr-type gravitomagnetic potential) decreases 
monotonically from (z —>■ 00) to some constant value (z — ► —00). We 
have assumed the direction of (3 along x-axis. This shift function derives an 
MHD wind which extracts translational energy analogous to the rotational 
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energy for the Kerr metric. The shift vector in three dimensions will be 
denoted by the Greek letter {3. The Kerr-type horizon has been pushed off 
to z — — oo. The pair production region lies at z = where the plasma 
is created. The newly created particles are then driven up to relativistic 
velocities by magnetic-gravitomagnetic coupling as they flow off to infinity 
and down towards the horizon. Geometrized units will be used throughout 
the paper. 



2.2 Electrodynamics in Kerr Planar Analogue 

We consider the magnetosphere filled with MHD fluid and take the perfect 
MHD flow condition in fluid's rest-frame 



E + VxB = (2.3) 

with V, B and E are fiducial observer (FIDO) measured fluid velocity, mag- 
netic and electric fields respectively. For perfect MHD flow in (12 .11) with 
a = 1, differential form of Faraday's law in 3+1 formalism [18] turn out to 
be 

+ (B.V)/3 - (V./3)B = Vx (VxB), (2.4) 

GST 

where = — /3.V is the FIDO measured rate of change of any three- 
dimensional vector in absolute space. Gauss law of magnetism according to 
FIDO can be written as [18] 

V.B = 0. (2.5) 

For (12.11) with a = 1, the local conservation law of rest-mass [18] according 
to FIDO is 

D,, ° + p 7 2 V.^ + Po V.(V - 0) = 0, (2.6) 



Dt ' ' Dt 

Dt ~ dr 



where po is the rest-mass density, 7 is the Lorentz factor and -j^- = -f- + 
V.V = + (V — /3).V is the time derivative moving along the fluid. The 
FIDO measured law of force balance equation [18] for the spacetime, given 
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by Eq. fl2.ip with a = 1, takes the form 




(2.7) 



where [l is the specific enthalpy and p is the pressure of the fluid. The FIDO 
measured local energy conservation law (Eq.(2.4) of |28j), for Eq. (12. II) with 
a = 1, is given as follows 



Eqs.(j23D-d2lD give the perfect GRMHD equations. 

3 Specialization of Background Flow for Model 
Spacetime 

In this section, we give the background flow and relative assumptions which 
will be used to simplify the problem. 

3.1 Description of Flow Quantities 

The FIDO measured 4-velocity of fluid can be described by a spatial vector 
field lying in the xz-plane [Tj5] 




+ Po /i7 2 (V.V) - Po /i7 2 V.(V.V)/3 + ^{(V x B).(V x B) 

+(V x B).-^-(V x B) + (V x B).{(V x B).V}/3 

(It 

-(V x B).(V x B)(V./3)} = 0. 



(2.8) 



V = V(z)e* + u(z)e z . 
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Here the Lorentz factor takes the form 7 = / 1 _^_ V 2 - The magnetic field 
measured by FIDO is also assumed to be in the ^-direction 

B = 5{A(z)e x + e z }, 

where B is constant. The corresponding Poynting vector becomes 

S = — (E x B). 

We have considered an example of stationary flow of an isothermal MHD fluid 
in our model spacetime (12.21) . These flows are used as stationary model mag- 
netospheres whose dynamical perturbations are to be studied. The plasma 
is moving in the xz-direction. The perturbed flow is along ^-direction due 
to the black hole's gravity and along x-direction due to rotation of the black 
hole (in direction of shift vector of our analogue spacetime). This flow will 
be analyzed to seek the properties of plasma waves. 

3.2 Perturbations and Wave Propagation 

The perturbed flow in the magnetosphere (which is in the xz-plane) can be 
characterized by velocity V, magnetic field B, the fluid density p and pressure 
p. We denote the unperturbed quantities by a superscript zero and the fol- 
lowing dimensionless notations are used for perturbations (<5V, 5B, 5p, 5p) 

dp „, 5p . 

P = — = p{t,x,z), p= — =p{t,x,z), 
P V 

v = 5V = v x (t,x, z)e x + v z (t,x,z)e z , 
b = = b x (t,x,z)e x + b z (t,x,z)e z . (3.1) 
The perturbed variables take the following form 

P = P° + PP, P = P° + m, 

V = V° + v, B = B° + 5b. (3.2) 

It is also assumed that the perturbations have sinusoidal dependence of t, x 
and z. Thus 

pit, X, Z) = Cie -i("t-k x x-k*z) ^ X; z j = c&e -i^t-k x x-k z z) ^ 

V x (t, X, Z) = C2e -i{^-k x x-k z z) ^ Xj ^ = C3e -i(ut-k x x-k z z) ? 

b x (t, x, z) = c i e- i( - ut - k ' x - k ' g \ b z (t, x, z) = C5e -i("t-k x x-k z z) _ 
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Using the values of components of k, we can discuss the quantities like phase 
velocity vector, group velocity vector, refractive index and its change with 
respect to angular frequency These quantities would help to investigate the 
wave behavior of the Kerr black hole magnetosphere and the properties of 
Veselago medium. 



4 GRMHD Equations for Kerr Spacetime in 
Isothermal State of Plasma 



The isothermal equation of state means that there is no exchange of energy 
between the plasma and the magnetic field. This state can be expressed by 
the following equation 

P + V 

\i = = constant. 

Po 



When we use this equation of state, the set of GRMHD Eqs. (l2l| -( |2T8]) take 
the following form for the spacetime given in Eq. fl2.2p . i.e, f3 = (/3 X , 0, 0) 



dB 



+ (B.V)/3 = Vx (VxB), 



dr 
V.B = 0, 

D(p + p) 



DY 



Dr +<P+P)1? V - : DT 



+ (p + p)V.V = 0, 



(4.1) 
(4.2) 

(4.3) 



{ (P + Ph 2 + *V + (P + ptfViVj - 1-BiBi 

Hp + v)i 2 Vi, k v k + ( P + vHViVjV\ k v k = ( P + 



dV j 
dr 



-P.< 



4tt I V dr 



(4.4) 



7 2 V.^(p + p) + 2(p + p) 7 4 V.^ - ^ + (p + p) 7 2 ( V.V) 
Dt Dt dr 

-(p + p) 7 2 V.(V.V)/3 + ]- [(V x B).(V x B) 

47T 



+(V x B).-^-(V x B) + (V x B).{(V x B).V}/3 

dr 



0. 



(4.5) 



These equations proceed in a similar way as used in [27j-[28j. Equations 
(13.11) and (13.21) as well as the restrictions for the velocity and magnetic fields, 
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given in Section 3.1, lead to the perturbed form of Eqs. (14 . 1 j) - ( 14 . 5j) given in 
Appendix A. When we use Eq. ( 13.31) . the Fourier analyzed perturbed equa- 
tions take the following form 

ikzCi — (tk z X + A')c3 — Ci{ik z u — loj + u) 

+c 5 {(V - (3)ik z + (V- (3)'} = 0, (4.6) 

ik x C2 — ik x \cz + lcs{(V — (3)k x + k z u — u} = 0, (4.7) 

k x c A = -k z c 5 , (4.8) 

d[Lp{-uj + (V - (3)k x + uk z } - {p'u + pu' + pwy 2 (W + uu')}} 

+c 2 (p + p) [-i(u + (3k x )-f 2 V + ik zl 2 uV + tk x (l + 7 V 2 ) 

+7 2 m{(1 + 2 1 2 V 2 )V + 2 1 2 uVu'}] + c 3 (p + p) [-l(u + /3k x )^ 2 u + ik xl 2 uV 



u 



+tk z (l + 7 V) - (1 - 2 7 V)(1 + 7V)- + 2 7 4 wW 



u 



+ C 6 {ip{-UJ 



+{V- p)k x + uk z } + {p'u + pu' + pu-f 2 (W + uu')}} = 0, 



(4.9) 



cip 7 2 u{(l + -? 2 V 2 )V + ~f 2 uVu'} + c 6 p[ 7 2 M{(l + j 2 V 2 )V + ~f 2 uVu'} + ik x \ 
+ (3k x ) ((p + p) 7 2 (l + 7 2 V 2 ) + B ~ 



+c 2 



47T 



B 2 



+i{k x V + k z u) { (p + p) 7 2 (l + 7 V 2 ) - — } + (p + p) 7 4 w{(l + 4rfV 2 )uu' 



+4(1 + 7 V 2 )VK'}] +c 3 
+i(k x V + k z u) \ (p + p)-f\V + 



-n,,j + (3k x ) { (p + p)- i aV - 
\B 2 



\B 2 



An 



4tt 

+ (p + p) 7 2 [2 7 2 (l + 2 7 V)wVV 



B 2 uX' 
An 



+{(1 + 2 7 V)(1 + 2YV 2 ) - YV 2 }V] + 



d2 n2 

+— cA-ikJl - it 2 ) + W} + — c 5 {-A' - u(V - P)' + ikJl - V 2 
Air An 1 

-2iuVk z } = 0, 

ci 7 2 p[m{(1 + 7 V)u + 7 2 lW} - V r /3 / ] 
+c 6 [ 7 2 p{(l + 7 V)m/ + 7 VwV - V(3'} +p + ptk z ] 

XB 2 



(4.10) 



+c 2 



t(o; + /3^) <^ (p + p) 7 W 



47T 



+ t(k x V + A; 2 m) { (p + p) 7 4 w1/ 



A£? 2 
An 



+ (p + p) 7 2 { 7 W(l + 4 7 V 2 ) - /?'(! + 2 7 V 2 ^ 



+2V 7 W(1 + 2 7 V)}] +c 3 [-L{uo + (3k x ) {(p + p) 7 2 (l + 7 V 

+ i{KV + k z u) j(p + p) 7 2 (l + 7V) - 
+ (p + P)7 2 K(1 + 7V) (1 + 4 7 V) + 2w 7 2 {(l + 2 1 2 u 2 )W - V(3'}} 



r>2 

•t-AA'u 

4tT 



i? 2 Z? 2 

+ — c 4 {iA; 2 A(l - it 2 ) + A' - AW} + — c 5 {2XuVik z 

4:71 4:71 



+Xu(V - (3)' - Xik x (l - V 2 )} = 0, (4.11) 
Ci 7 2 [— mop + p'u + pu + 2pw 7 2 (l / l // + uu) — fq 2 uV(3' + ik x p{V — (3) 
+pik z u] + c 6 [— iup(-y 2 — 1) + 7 2 {p'u + pu' + 2pw 7 2 (l / l // + m/) 
-P7 2 m1//?' + - (3) + pik z u}} + c 2 [-^{2(p + p) 7 V 

-^(«A - V)} + ik x [(p + p) 7 2 {l + 2 1 2 V(V - f3)} - ^(V - 0){uX - V)} 

B 2 

+tk z u{2{p + p) 7 4 V + — (uX - V)} + (p + p) 7 \{2 7 V 

47T 

+6 7 V(W' + uu') - (3'{1 + 2 7 V 2 )} - -— ] + c 3 [-iu{2(p + p) 7 4 M 

47T 

A52 («A - V)} + tk x {V - (3){2(p + pHu + ^{uX - V)} 



4tx " v 4vr 

+^{(p + p) 7 2 (l + 2 7 V) - ^{uX - V)} + (p + p) 7 2 {-- 

47T M 

+2 7 2 W + 67V (W + uu) + ^{W + uu) - V(3'{1 + 7 V)} 

{A - u(uX - V)}} + c 4 —[u{X' - (uX - V)u'} 
4tt 4tx 

B 2 

+ L k z {uX -V)(l- u 2 )} + c 5 — {-X'V + u(uX -V)(V- (3)' 

4n 

-ik x (uX - V)(l - V 2 ) + 2tk z uV} = 0. (4.12) 

Eq. fl4.8p gives the relation between x and ^-components of the wave vector 
i.e., k 7 = ——k T which will be used in the next section. 



5 Numerical Solutions 

This section is devoted to the numerical solutions of the dispersion equations. 
The following subsection contains the relative assumptions which make the 
equations easier to deal with. 
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5.1 Relative Assumptions 



In our stationary symmetric background, the x-component of the velocity 
vector can be written in the form [19J V = C + Ait, where C = (3 + Vp with 
Vp as an integration constant. We assume the value of the shift function [19] 
(3 = tanh(z) — 1 with Vp = 1. Further, B 2 = 8ir and A = 1 are taken so that 
the magnetic field becomes constant. Thus the x-component takes the form 



Substituting the value of V in the conservation law of rest-mass for three- 
dimensional hyper surf ace, i.e., 



with the assumption that rest-mass density is constant, we get an equation 
of the form 

3u 2 + 2utsmh(z) + tanh 2 (» -1 = 

quadratic in u with the assumption that A/p = 1. The solutions of this 
equation and the corresponding values of V are given as follows 



We shall use these values to solve the dispersion relations. The Poynting 
vector for these values takes the following form 



These quantities are valid for the region outside the pair production region. 
We consider the region —5 < z < 5 and omit the region — 1 < z < 1 due to 
large variations in the background flow quantities. In the rest of the region, 
these quantities become constant and Fourier analyzed procedure is valid for 
this region. Further, we use the relation k z = —k x which reduces the wave 
vector to (k x , 0, —k x ). 



V = 1 + p + u = tanh(z) + u. 



Po7« = Mp + v)l u — A (constant) 



(5.1) 




(5.2) 



(5.3) 



S = 2 tanh(z) (e x — e z ) . 
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The computer programming (using Mathematica) is used to evaluate a 
root of the dispersion relation for the plasma moving towards the black hole 
with the velocity components given by Eq. fl5.2p . This is given as a separate 
file with all the required codes. Other roots can be obtained in a similar 
manner. 

It is observed that the sextic equation has all roots admitting imaginary 
values at several points. The quintic equation gives one real root of the 
positive z region for both the values of velocity (Eqs. (15.21) and (I5.3P ) shown 
in the Figures 1 and 2. For the negative z region, the velocity, given by 
Eq. (l5.2p . gives one real root shown in the Figure 3. The approximated root 
becomes imaginary at z = —5 which we omit and our mesh reduces to —4.8 < 
z < —1, < uj < 10 for the interpolating function. The values at z — —5 
are extrapolated afterwards. The velocity components, given by Eq. (l5.3p for 
negative z region, leads to three real roots shown in the Figures 4, 5 and 6. 
The real data values for the root give a real interpolation function. 

It is clear that the Figures 1 and 2 represent the neighborhood of the pair 
production region towards the outer end (as the wave number is found for the 
positive values of z) whereas the Figures 3, 4, 5 and 6 show the neighborhood 
of the pair production region towards the event horizon (as the wave number 
is determined for the negative values of z) . 

5.2 Results 

First, we obtain k x for the velocity components, given by Eqs. (15.21) and (15.31) 
in the positive z region. These are shown in the Figures 1 and 2 respectively. 

In the Figure 1, the x-component of the wave vector is negative near the 
pair production region and for the waves with negligible angular frequency. 
For each angular frequency, the waves grow monotonically in number when 
they move away from the event horizon. There is a sudden increase in the 
x-component of the wave vector, it admits positive values for a particular 
value of z, then decreases a bit and smoothly increases afterwards. The fluid 
near the pair production region (region with negative x-component of the 
wave vector) possesses the negative values for the x-components of the phase 
and group velocities. For this region, the wave vector is in the opposite 
direction of the Poynting vector and hence it shows the existence of Veselago 
medium there [23]. For the same region, the phase and group velocity vectors 
are in the direction opposite to the Poynting flux and hence the regions 
are of negative phase and group velocity propagation. The change in the 



12 



x-component of Wave Vector 




Figure 1: For the velocity components given by Eq. ( 15.21) . plasma admits 
the properties of Veselago medium near the pair production region. As the 
waves move away from the pair production region, they disperse normally. 
Negative phase and group velocity propagation regions are observed near the 
pair production region. 



13 



refractive index with respect to the angular frequency is positive for the 
regions (i) 1.6 < z < 2, 0.4 < u < 1, (ii) 2 < z < 3, 0.42 < u < 2.92, (iii) 
3 < z < 4, 0.225 < w < 10 and (iv) 4 < z < 5, 0.14 < uj < 10 for which the 
dispersion is normal [26], [30]. In the rest of the region, most of the points 
admit anomalous dispersion. 

The Figure 2 shows that the x-component of the wave vector is negative 
for the region 1.0 < z < 1.89. It is large near the event horizon and decreases 
up to z = 1.75 after which it increases and fluctuations occur in the values. 
In the region 1.89 < z < 10, it takes random values. The negative values of 
k x in the region implies that the wave vector is in the opposite direction to 
the Poynting vector which indicates the properties of Veselago medium. In 
the region < z < 1.89, the rc-components of the phase and group velocities 
take negative values and hence this is of negative phase and group velocity 
propagation region. Both these quantities admit random values in the region 
1.89 < z < 10. For the region 1 < z < 1.6, < uj < 0.079, the change in the 
refractive index with respect to the angular frequency is positive and hence 
the dispersion is found to be normal. In the region 1 < z < 1.6, 0.079 < 
uj < 10, the quantity ^ < which indicates anomalous dispersion in this 
region [26]. The rest of the region shows random points of normal as well as 
anomalous dispersion. 

For the negative z region, i.e., the region towards the event horizon of the 
black hole in the neighborhood of the pair production region, we obtain one 
value of k x for the velocity components given by Eq. fl5.2p and three for the 
velocity components given by Eq. (15.31) . These values are shown respectively 
by the Figures 3, 4, 5 and 6. 

In the Figure 3, the x-component of the wave vector is negative for the 
region —1.925 < z < —1.0 where the Poynting vector is parallel to the 
wave vector and hence the medium is usual. The refractive index greater 
than one and positive change in the refractive index with respect to the 
angular frequency indicate normal dispersion. In the rest of the region, all 
the three quantities admit random values and hence there are normal as well 
as anomalous points of dispersion. For the region —1.4 < z < —1.0, 0.5 < 
uj < 10, the x-components of phase and group velocities are negative such 
that v px > v gx . These velocity components admit random values in rest of 
the region. 

The Figure 4 indicates that the x-component of the wave vector is nega- 
tive for the whole region and hence the wave vector and the Poynting vector 
are in the same direction showing the existence of the usual medium. As the 
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x-component of Wave Vector 




Figure 2: Veselago medium exists near the pair production region for the 
velocity components given by Eq. fl5.3p . In the same region, negative phase 
and group velocity propagation regions are observed. Most of the region in 
the neighborhood of the pair production region shows anomalous dispersion. 
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x-component of Wave Vector 






Figure 3: Near the pair production region, the dispersion is normal whereas 
the rest of the region admits normal as well as anomalous points of dispersion 
for velocity components given in Eq. (15.21) . 
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x-component of Wave Vector 




values of z and uj grow, k x decreases and hence v px and v gx are negative in 
this region. Although v px > v gx for the region < uj < 10~ 15 , yet they are 
nearly equal for rest of the region. The refractive index is greater than one 
in the whole region. The refraction increases as the waves move towards the 
pair production region. The change in the refractive index with respect to 
the angular frequency is negative for the region < uj < 0.6 which shows 
that anomalous dispersion of waves. For the rest of the region, it is positive 
and thus indicates normal dispersion. 

The Figure 5 shows that k x is positive for the whole region. Thus the wave 
vector is in the opposite direction to the Poynting vector which indicates the 
presence of Veselago medium. The quantity k x increases with the increase 
in z and uj except for the waves with negligible angular frequency. v px and 
v gx are nearly equal and admit positive values which show negative phase 
and group velocity propagation in the whole region. The refractive index 
is negative and decreases as the values of z and uj increase. The change in 
the refractive index with respect to the angular frequency is negative for the 
regions (i) -2 < z < -1, 2.15 < uj < 10 (ii) -3 < z < -2, 5.5 < uj < 10 and 
(iii) —4 < z < —3, 9.5 < uj < 10 which indicates anomalous dispersion in 
these regions. It is positive for the rest of the region which indicates normal 
dispersion except for the waves with negligible angular frequency. 

In the Figure 6, the x-component of the wave vector is positive for the 
whole region and increases with the increase in the angular frequency. As 
the waves move away from the pair production region, the re-component of 
the phase velocity decreases slightly and then increases. In contrast, the x- 
component of the group velocity increases a little and then decreases. The 
refractive index is negative due to the fact that the Poynting vector is in the 
opposite direction to the wave vector which shows the existence of Veselago 
medium. The positive values of x-components of phase and group velocities 
show the existence of negative phase and group velocity propagation regions. 
The change in the refractive index with respect to the angular frequency is 
negative throughout the region and hence shows anomalous dispersion except 
for the waves with negligible angular frequency. 

6 Conclusion 

It is well-known that charged particles are created in the pair production 
region. Some of these particles which are pushed on to orbits with negative 
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x-component of Wave Vector 




x-component of Phase Velocity Vector x-component of Group Velocity Vector 




Figure 5: For the velocity components given by Eq.f l5.3p . Veselago medium 
exists in the whole region with negative phase and group velocity propagation 
property. The region of normal dispersion extends as the waves move away 
from the pair production region. The waves with negligible angular frequency 
do not disperse normally. 
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x- component of Wave Vector 




Figure 6: For the velocity components given by Eq. (15.31) . plasma admits the 
properties of Veselago medium. The dispersion is anomalous except for the 
waves admitting negligible angular frequency. The negative phase and group 
velocity propagation regions are also observed. 
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energy by the Lorentz force move towards the event horizon and others move 
towards the outer end of the magnetosphere. These particles would reach 
their destinations if the plasma existing in the neighborhood of the pair pro- 
duction region allows them to do so. The generation of plasma is necessary to 
support the MHD magnetically dominated flow. Due to particle generation, 
waves are produced in the neighborhood of the pair production region. The 
dispersion relations of waves lead to understand how much the surrounding 
medium let the waves to disperse through. 

This paper studies the wave properties for the isothermal plasma moving 
with velocity V and admits a constant magnetic field which thread the Kerr 
black hole magnetosphere. The gravitomagnetic waves and the pair of parti- 
cles are produced in the z = region. If the medium living around the pair 
production region allows the particles and waves to pass through, the energy 
extraction from the black hole is possible. This can be well understood by 
investigating properties of the waves in this region. 

We have considered a black hole immersed in a rarefied plasma with 
uniform magnetic field which seems to provide support for carrying currents 
flowing across the magnetic field lines. Due to the strength of the magnetic 
field, a lot of energy can be extracted due to the plasma particles that fall into 
the black hole's horizon has negative energy. The 3+1 GRMHD equations are 
derived for the neighborhood of pair production region and two-dimensional 
perturbations are discussed in the context of perfect MHD condition. We 
assume that the rotation is in the x-direction and the horizon is at z = 
— oo. The perturbations are taken only in the ^-direction. The dispersion 
relations are formulated by assuming the perturbations as plane waves. We 
solve these relations by taking the wave vector as (k x , 0, — k x ) and obtain the 
re-component of the wave vector. This component leads to properties of the 
isothermal plasma in the neighborhood of the pair production region. 

We have discussed these relations for the regions 1 < z < 5 and —5 < 
z < —1. The dispersion relations for the region 1 < z < 5 are shown in 
the Figures 1 and 2. These Figures indicate that near the pair production 
region, the plasma admits the properties of Veselago medium. The region 
which is nearer to the pair production region does not allow the waves to 
pass through. Thus the particles and waves cannot get out of this region. 
The small regions far away from the pair production region admit normal 
dispersion of waves which indicate that the waves pass through them. As we 
go away from the pair production region, normal dispersion exists frequently 
as shown in the Figure 1. 
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The region —5 < z < — 1 allows us to investigate whether there is a 
possibility for the waves to move towards the black hole event horizon or 
not. The dispersion relations for this region are shown in the Figures 3, 4, 5 
and 6. From the Figures 3, 4 and 5, we find that there are chances for the 
waves to pass through the neighborhood of the pair production region when 
the plasma admits the properties of usual or Veselago medium. The Figure 
6 indicates that there can be situation when plasma admits the properties of 
Veselago medium in the neighborhood of the pair production region, it may 
not allow the waves to pass through the region. 

It is interesting to note that the Figures 2 and 3 show the irregular depen- 
dence of wave vectors on the angular frequency and z. Mathematically, this 
irregularity is due to the nature of the roots obtained for these graphs. The 
irregular behavior may be due to a disturbance of the equilibrium between 
outward and inward directed forces. The outward directed forces are caused 
by the particle pressure and the curvature drift due to non-uniform magnetic 
field and inward directed forces are exerted by the tangential stress of the 
magnetic field lines for the low frequency regime. 

For the high frequency regime, there is a class of MHD instabilities which 
sometimes develop in a thin plasma column carrying a strong current. If a 
kink begins to develop in such a column, the magnetic forces on the inside of 
the kink become larger than those on the outside, which leads to the growth 
of perturbation. The column then becomes unstable and causes a disruption. 
Both the ballooning and kink modes are ideal MHD instabilities. 

In the Figures 1, 2, 5 and 6, we obtain the properties of Veselago medium. 
The phase and group velocity vectors propagate in the direction opposite to 
the Poynting vector which verify the results of Mackay et al. [22] according 
to which rotation of a black hole is required for the negative phase velocity 
propagation. 

We can conclude that waves produced in the pair production region due 
to pair creation cannot get out of its neighborhood towards the outer end of 
the magnetosphere. The same result has been obtained for the cold plasma 
case [22] • We obtain some cases where favorable conditions are present to 
allow energy to move towards the black hole horizon. For the cold plasma, 
these conditions are present for the usual medium whereas for the isothermal 
plasma, these situations occur for usual as well as Veselago medium. For the 
plasmas with pressure, the black hole can suck particles and waves for both 
the usual and Veselago medium whereas for the cold plasmas, this situation 
holds for the usual medium. 
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The strong magnetic coupling enforce the accreting particles to fall into 
the black hole with negative energy and negative angular momentum. This 
indicates that energy and angular momentum flow from the black hole into 
the disk. When the particles fall into the black hole with negative energy, 
energy is extracted from the black hole [31]. When the particle with positive 
energy and positive angular momentum leaves the pair production region 
and goes towards the event horizon, much energy and momentum are lost 
and ultimately the particle has negative energy and angular momentum [32J . 
Thus if the particle either with negative or positive energy leaves the pair 
production region and gets a chance to reach the event horizon, the result 
is the extraction of energy from the black hole transmitted to the accretion 
disk. This shows that when the magnetosphere is filled isothermal plasma 
admitting the properties of Veselago as well as usual medium, our results 
indicate that energy extraction is possible. 
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Appendix A 

This Appendix includes the details to reach at the perturbed form of the 
GRMHD equations (14. ip ~ (14.51) . The component form of these equations is 
also given. 

When we introduce the perturbations from Eq. (I3.2j) . the linearized GRMHD 
Eqs.flO)- (l43|) become 

d 1 

— - p.V){5B) >=Vx(vxB) + Vx(Vx <5B) — (<JB.V)/3(A1) 



dt 

V.(5B) = 0, (A2) 

+ ( V - 0).v} (5p + 6p) + (p + ph 2 V. | ^ + (V - /?). v} v 

-(p + p)(V.v) + (dp + 5p)(V.V) + (Sp + 5p) 7 2 V.(V.V)V 
= -2(p + p) 7 2 (V.v)(V.V) In 7 - (p + p) 7 2 (V.VV).v 
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+ (p + p)v.Vlnw, (A3) 



Up + Ph 2 + ^| hi + <J> + pft'ViVj - ^B,B, 



at "- v ^ 



Bx(Vx«l 



dr 



Hp + ph\jV j + (p + P h 4 v t v hk vw k + i- 

-i- {(5^), - # = -(6p) ti + i 2 [(5p + <fp) V' 

+2(p + p) 7 2 (V.v)^' + (p + pK]/^ + - fl,- 0*5' 

-(p + p) T 4 (^ + ^)l4 ) .^ 

- 7 2 {(<5p + 5p)V J + 2(p + p) 7 2 (V.v)V^ + (p + p)^'}^- 

-7 4 ^{(5p + 5p)V j + 4(p + p) 7 2 (V.v)\/ J + (p + p)^'}^, (A4) 

7 2 + (V - /3).v| (5p + <5p) + 2(p + p) 7 4 V. j-| + (V - /3).v| v 

-(P + p)7 2 v.V In u + (p + p) 7 4 V.(v.V)V + 6(p + p) 7 6 (V.v)V.(V.V)V 
+2(5p + 5p) 7 4 V.(V.V)V + 2(p + p) 7 4 v.(V.V)V 

-2(p + p) 7 4 (V.v)V. V In u + 2(p + p) 7 4 (V.v)(V.V) - ^(5p) 

+ (5p + 5p) 7 2 (V.V) + (p + p) 7 2 (V.v) - 7 2 03.V)(<5p + 5p) 
+2(p + p) 7 4 (V.v)03.V In u) - 6(p + p) 7 6 (V.v)V.(/5.V)V 
-2(p + p) 7 4 v.(/3.V)V - 2(5p + 5p) 7 4 V.(/3.V)V 
-(5p + 5p) 7 2 V.(V.V)/3 - (p + p) 7 2 V.(v.V)/3 

-2(p + p) 7 4 (V.v)V.(V.V)/3 + ^[(vxB).(Vx B) 

47T 

+(V x 5B).(V x B) + (V x B).(V x SB) 
+(V xB).{-xB + Vx— j 



0. (A5) 



The component form of these equations are 

^7 + Vb X:X + = -u'b x + (V - (3)'b z + v x , z - Aw 2jZ - A'u z , (A6) 

^7 + V& Z)a . + ub Z:Z = \v Z:X - v X:X , (A7) 
b x , x + b z , z = 0, (A8) 
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dp dp T _ T _ _ _ _ _ N . . 

p— + p— + pi/ p >a; + pv/p,^ + pwp iZ + (p p) \P u ~\~ pu 
dr dr 

+pu 1 2 {W + uu')} + (p + p) 7 2 (V ^ + + (p + p)(l + 7 V 2 )^ 

+ (p + p)(l + 7 2 u 2 )v z , z + (p + p)u^7 2 (u x ,« + v z , x ) 
= -(P + pW u \( 1 + rfV 2 )V' + 2 1 2 uVu'\v x 



If 

+ (p + p)[{\ - 2 7 V)(1 + 7 V)- - 2 1 4 u 2 VV']v z , (A9) 

+ | (p + p) 7 2 (l + 7V 2 ) - — I (V^ + ut;^) - -^uVb gtZ 

+ |(p + p)7 4 «V + (^ + uv z , z ) + ^{(1 - V 2 )^ 

Z? 2 i? 2 

-(1 - u 2 )b x , z } = -^uu'b x + —{A' + u(V - /3)'}6 a - pp, x 

~(PP + PPhM(! + 1 2 V 2 )V + 7 2 wW} - (p + p)7 4 w{(l + 4 7 V 2 )W 
+4(1 + *fV 2 )VV'}v x - [(p + p)7 2 [{(l + 2 7 V)(1 + 2 7 V 2 ) - -f 2 V 2 }V 

B 2 

+2 7 2 (1 + 2 7 V)wlV] + — u\']v z , (A10) 

47T 

{ (p+ph2(1+7 v ) + ^!}^ + {( P+ ^-^}f 

f X 2 B 2 1 A5 2 

+ Up + p)7 2 (l + 7 2 ^ 2 ) - I (Vu*,* + uv z , z ) + —^^-uVb ZjZ 

( \B 2 1 A-B 2 

+ I (p + p) 7 4 ^ + \ (Vv XtX + uv X!Z ) - ^{(1 - V 2 )b z>x 

B 2 \B 2 
-(1 - v?)b x , z } = -—(A' - AW)fo s - -;— «(V - 

-PP7 2 {W(1 + 7V) + j*u 2 VV' - V(3'} - [pp i2 + p'p 

+pp 1 2 {uu\l + 7 V) + 7 V VV - V{3'}} - (p + p) 7 2 [ 7 V(l + ^ 2 V 2 )V' 

-(1 + 27 2 \/ 2 )/5' + 2 1 2 uV{\ + 2 7 V)wV - [(p + p)7 2 {-27 2 M\//5' 

+ (1 + 7 2 m 2 )(1 + 4 7 Vy + 2 7 2 (1 + 2 1 2 u 2 )uVV'} - —— u\']v z , (All) 



47T 



7 2 P? + P(7 2 - 1)^ + P7 2 {A + pu' + 2pu~i 2 {VV + uu') 
at at 
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-puV(3'} + p7 2 { up' + u'p + 2pwy 2 {W + uu') - p-f 2 uVf3'} 
+j 2 pp, x {V - 0) + l 2 upp :Z + i 2 pp, x {V - (3) + i 2 upp, z 

f)u flu \ 

+ ^{2(p + p) 7 V - — (uA - 1/)} + -^{2{p + p) 7 4 n + —{u\ - V)} 

+v x>x [(p + p) 7 2 {l + 2 T V(y - (3)} - ^{V - (3){uX - V)] 

+v x>z u[2(p + Ph 4 V + ^{uX- V)} + v ZiX {V - (3) [2(p + p) 7 4 « 

A7T 

+^( uA _ V)] + „ { ( p + p )(i + 2 7 V) - ^u(uX - V)} 

47T " 47T 

+^_( mA - - u 2 )b x , z - (1 - F 2 )^,, + 2uFfo, )Z )} 

ATT 

+— w6 x {A' - (uX - V)u'\ + — b z {-X'V + u(uX -V)(V- (3)'} 

ATT All 

B 2 X' 

+v x [{p + p)i 2 u{2 1 2 V + 6 7 V(\/V r/ + uu') - 13' {1 + 2 7 V 2 )} - — - ] 

An 

+v z [— {X - u(uX - V)} + (p + ph 2 {-- + 2 7 W + 6 7 V( W + uu') 
An u 

W(VV + uu) - V(3'{1 + 7V)}] = 0. (A12) 

We have used the conservation law of rest-mass for three-dimensional hyper- 
surface, i.e., given by Eq. fl5.13p to simplify Eq. (IA9p . The same equation will 
be used to simplify the Fourier analyzed form of Eqs (4. 11)- (4. 17). 
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